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Abstract 

This paper is a contribution to the Proceedings of the Workshop 
Complexity, Metastability and Nonextensivity held in Erice 20-26 July 
2004, to be published by World Scientific. We propose a generaliza- 
tion to Merton's model for evaluating credit spreads. In his original 
work, a company's assets were assumed to follow a log-normal pro- 
cess. We introduce fat tails and skew into this model, along the same 
lines as in the option pricing model of Borland and Bouchaud (2004, 
Quantitative Finance 4) and illustrate the effects of each component. 
Preliminary empirical results indicate that this model fits well to em- 
pirically observed credit spreads with a parameterization that also 
matched observed stock return distributions and option prices. 

1 Introduction 

Recent years have brought with them a vastly growing interest in 
credit markets and the related credit derivatives markets. The amounts 
trading on these markets are growing rapidly. As the liquidity has in- 
creased, so has the interest in models which can price credit risk. 



1 



Indeed, one of the most esteemed pricing models was proposed al- 
ready in 1974 by Merton [1], although it is only recently that the 
credit markets have become liquid enough and transparent enough to 
actually use these models to price the yield on risky corporate bonds, 
whose spread to Treasuries is also called the credit spread. It is quite 
interesting to note, as Miller and Modigliani pointed out in 1958 [2], 
that there is a relationship between a firm's debt and its equity in 
that the sum of the two arc equal to the total assets of the com- 
pany. Therefore, a company should be indifferent to how it finances 
its projects, via debt or equity. On the other hand, the equity mar- 
kets and bond markets are so much more developed than the credit 
and credit derivatives markets. This discrepancy together with the 
intimate relationship between the different markets implies that there 
could be arbitrage opportunities if relative mispricings exist between 
the fair value of the debt relative to the fair value of the equity. Such 
a capital structure arbitrage has been a very popular trading strategy, 
especially in the early 2000's. 

In reality, the capital structure of a company is quite complex. For 
example, corporate debt has different levels of seniority, which defines 
the order in which the debt holders get paid. Clearly, corporate debt is 
risky since there is a chance that the company will default before being 
able to pay back its debt. Equity holders get paid after all debt has 
been repaid, but even among stock holders there are different priority 
levels. However, as Merton pointed out [1], in the simple case where a 
company has just one type of debt outstanding, the stock price can be 
seen as a call option on the total underlying assets, struck at the face 
value of the debt. Likewise, the corporate bond can be interpreted 
as long a riskless bound plus a short put on the assets, struck at the 
face value of the debt. Therefore, using the Black, Scholes and Mer- 
ton option pricing theory [3, 4], it is possible to price corporate debt 
with the techniques of option pricing. Merton's structural model, and 
similar models which derive from Morton's ideas, constitute a widely 
used framework for pricing debt and credit risk, and for obtaining 
estimates of the risk-neutral probability of default of a company. 

The basic strength of Merton's model lies in the ability to price the 
debt as an option on the underlying asset process, via standard tools 
of option pricing theory. The approach is complicated by the fact that 
in reality, the asset process is not observable and must be backed out 
from the stock process. It is typical that a log-normal distribution is 
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assumed for the asset returns, and that the driving noise for both the 
asset and the stock dynamics follows a standard Brownian process. 

However, it is quite well-known that the empirical distribution of 
log returns is highly non-Gaussian, exhibiting fat tails which are ne- 
glected in the standard Black, Scholes and Merton option pricing the- 
ory. Indeed, in options markets one observes that the standard theory 
underestimates the prices of options at strikes below and above the 
current stock prices. This means that one must use a higher volatility 
parameter in conjunction with the Black-Scholes- Merton theory in or- 
der to correctly price the options. A plot of this volatility versus the 
strike price generally forms a concave function, rather than a straight 
line which would be the case if the model was perfect. There have been 
several attempts in the literature to accommodate this fact, including 
stochastic volatility models, multifractal models, local volatility mod- 
els, models that assume fat-tailed random noise such as Levy noise, 
GARCH-like models and recent multi-timescale models (for an up-to- 
date review see [5]). However, these models are often quite complex 
and the simplicity of the Black-Scholes- Merton approach is lost. A 
unique martingale measure is typically not found, neither are closed- 
form solutions. 

An alternative approach which one of us proposed recently relics 
on modeling the noise as a statistical feedback process [6, 7, 8] which 
yields non-Gaussian distributions for the stock returns yet maintains 
many useful features of the standard theory. In particular, closed form 
solutions which generalize the Black-Scholes-Merton model are found. 
These incorporate both fat tails and skew, two features of real returns 
which are absent in the standard theory. 

Since the non-Gaussian theory has shown some success in being 
able to price options on the underlying equity in a parsimonious fash- 
ion which matches well to empirical observations, our goal in this 
paper is to explore whether the same can be said about pricing credit. 
The basic notion is to extend Merton's model for pricing credit risk 
into the non-Gaussian framework. We shall then explore the effects of 
introducing fat tails and skew into the model. Finally we shall report 
some preliminary empirical results which indeed lend support to our 
approach. 
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2 Merton's model 



We proceed with a brief review of Merton's model. He developed a 
structural model relating the equity and debt markets by thinking 
of both stocks and bonds as being contingent claims on the same 
underlying, namely the assets of the company. His model assumed 
the assets A follow a standard log-normal process, 

dA = fiAdt + UAAdbj (1) 

where a; is a standard Brownian noise, b- correlated in time. As follows 
from the theorem of Miller and Modigliani, the value of the firm in 
invariant to its capital structure or in other words, that A is equal 
to the sum of its debt D and stock (equity) 5, such that the only 
relevant variable is the leverage L = D/A . Furthermore, one assumes 
that it is possible to continuously trade the assets. (Note that this is 
unrealistic in reality the stock and the bonds are traded, but not the 
assets themselves.) The interest rate r is assumed constant over time. 

In Merton's model, the scenario is such that the company has 
issued bonds of face-value D that will become due at a future time T. 
These bonds are assumed to be zero-coupon which means that there 
are no intermediate payments before expiry. In addition, the company 
has issued equity with no dividends payments. If at time T, the value 
of the assets A is less than the promised debt payment D, then the 
company defaults. The debt holders receive A < D and the equity 
holders receive nothing. If instead A > D, the debt holders get D and 
the stock holders get A — D. The equity of the company can therefore 
be seen as a European call option on the assets with maturity T and 
a strike price equal to D. In the standard framework this is just given 
by the closed-form Black-Scholes-Merton call option pricing formula 
(cf [9]), namely 

So = AoN{di) - De-'^N{d2) (2) 

where 

d, = i^^i^^^ + 0.5..VT (3) 

d2 = di-O.baAVT (4) 

and N{d) is the cumulative normal distribution calculated up to the 
limit d. Note that D = De~^^ is the present value of the promised 



4 



debt payment, so this expression can be expressed in terms of the 
leverage L = D/Aq, namely 

So = Ao{N{di) - LN{d2)) (5) 

This expression depends on Aq and the asset volatility a a, both of 

which arc unobscrvable. However, as shown by Jones et al (cf [10]), 
Ito's lemma can be used to determine the instantaneous asset volatility 
from the equity volatility leading to 

Sqcts^Aqua (6) 

Equations (2) and (6) allow the asset value Aq and the asset volatil- 
ity a A to be backed out as functions of the quantities So,as,T and 
L. Utilizing A = D + S, the present value of the debt can be thus 
calculated straightforwardly as 

Do = Ao- So (7) 

where Sq is given by Eq(2). Now, not that Dq can alternatively be 
expressed as 

Do = e-y^D (8) 

Namely as the promised debt payment discounted by some yield y. 
The difference between this yield y and the risk-free rate r defines the 
credit spread s on the debt, namely 

s = y-r (9) 

This is a very interesting quantity, because it can be shown to be 
roughly equivalent to premium on a the highly traded credit derivative 
called a Credit Default Swap (CDS). This is a transaction in which one 
party (the protection buyer) pays an annual premium to another party 
(the protection seller) who basically offers insurance on the bond of 
a company. If a default event (which is carefully defined) occurs, the 
protection seller delivers par to the protection buyer and receives the 
defaulted bond in exchange. (There are well defined criteria associated 
with this, somewhat similar to the delivery options in a bond futures 
contract. Cash settled versions of this contract also exist.) 
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3 A non- Gaussian approach 



The driving noise of the asset process in Merton's model is Gaussian, 
and the related equity process is also a Gaussian process with a con- 
stant volatility as which can be related to the asset volatility a a via Eq 
(6). Thus, the process for stock returns is a standard lognormal one, 
implying that the stock log-returns arc normally distributed across 
all time scales. Although this is a standard assumption in much of 
mathematical finance, leading to many interesting and useful results, 
it is surprisingly far from what one observes empirically. Indeed, real 
stock returns exhibit fat tails and peaked centers, only slowly converg- 
ing to a Gaussian distribution as the time scale increases [11, 12]. One 
model which generalizes the standard model in a way more consistent 
with empirical observations has recently been proposed by one of us 
[6, 7, 8]. In the spirit of that model, we generalize the asset price 
dynamics to follow a non-Gaussian statistical feedback process with 
skew namely 

dA = jXAdt + GAA'^d^l (10) 

dn = p^{n)dLo (11) 

Here a; is a standard Brownian noise, ij,a is the rate of return of the 
firm and cr is a variance parameter. The parameter a introduces 

an asymmetry into the model. 0, evolves according to a statistical 
feedback process [13], where the probability distribution P evolves 
according to a nonlinear Fokkcr-Planck equation 

dt 2 do? ' ^ ' 

This diffusion equation maximizes the Tsallis entropy of index q [14, 
15]. Equation (12) can be solved exactly, leading, when the initial 
condition on P is a P{rt,t = 0) = 5{Q), to a Tsallis distribution 
(equivalent to a Student distribution [16]): 

^=W)^^^^'^~ l)/3(t)J^^(i))"^ (13) 

with 

/3(i) = cH((2-g)(3-g)t)-^ (14) 

and 

Z{t) = {{2-q){3-q)cgt)^^ (15) 
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where the g-dependent constant Cq = -^{T'^i-^ — ^)) / {^'^ i-^)) ■ 

Eq. (13) recovers a Gaussian in the hmit g — 1 while exhibiting 
power law tails for all g > 1. The index q controls the feedback 
into the system. For q > I rare events (small P) will give rise to 
large fluctuations, whereas more common events (larger P) yield more 
moderate fluctuations. 

Pricing options based on this type of model was solved in previous 
work [7, 8]. For q = a = I, the standard Black-Scholcs-Mcrton model 
is recovered. For q = 1 but general a < 1, the model reduces to the 
Constant Elasticity of Variance (CEV) model of Cox and Ross (cf [8]). 
Closed form solutions for European call options based on this model 
were found. Evaluating the equity as a call option on the asset process 
struck at the debt D, we thus obtain 

So = Ao {1 + {1 - a)xiT)) — Pq{Qf)dnf 

J di 

- e'^'^D f^^ Pq{rtf)dnf (16) 

with Xf, d\ and ^2 defined in the Appendix. 

Just as in the standard Merton framework, it is possible to show 
that the relationship Eq (6) still holds in the non-Gaussian case. 

Though slightly more complicated, it is still possible to back out 
and a A from Eq (16) and Eq (6). The debt can then be calculated as 
Do = ^0 — 'S'o) with 5*0 as in Eq (16), and the credit spread can also 
be found as in Eq (9), namely 

Sq,a = y - r = -^log{^^) (17) 

The main difference is that now the spread is parameterized by q and 
a. Because Dq is based on So of Eq (16) allowing us to explore the 
effects of fat tails (^ > 1) and skew (a < 1). 

4 Results and Discussion 

Based on the standard Merton model, many practitioners would calcu- 
late the fair price of the credit spread from Eq (9). In order to obtain 
a value of s that matched the market value of the traded CDS, they 
would have to use a particular value of as- This implied as should in 
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Figure 1: Implied volatility of the credit spread as a function oi d = D/Aq for various 
values of q (the tail) and a (the skew). 

theory be the real volatility of the stock. In addition, it should the- 
oretically also be the correct volatility to price options of the stock. 
But in practice, it is already known that it is very difficult to find a 
single value of as that can price all stock options, let alone both stock 
options and credit spreads. 

In our previous work, we showed that indeed by introducing tails 
and skew into the underlying stock process, it was in many cases possi- 
ble to find a single value of as that matched well the underlying stock 
distribution, as well as observed stock option prices. The question we 
now ask is whether the values of q and a that fit well to the options 
market could also well-model observed spreads, or CDS prices. To 
answer this question we cite preliminary results [17] where the param- 
eters which best fit market spreads across a sample of 54 companies 
were implied based on a least squares fit between the theoretical spread 
values and market values. In that analysis, all the observables which 
enter the spread calculation were taken from reliable data sources as 
described in [17]. Certain assumptions were made in order to map 
the complex capital structure of real companies onto the simplified 
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single-debt scenario discussed here. The way in which this was done 
followed the lines of [10], where a weighted average of all outstanding 
debt was used to create a synthetic bond with a single face value and 
expiry. Interestingly, the parameter q was found to be in the range 
q = 1.2 to q = 1.4, and on average the parameter was a = 0.3 with lit- 
tle variation. These results are quite encouraging, since typical values 
that are found to fit option markets are very close to these values. An 
example is given in [8] where stock options on MSFT stocks are found 
to be well fit by g ~ 1.4 and a ~ 0.3, with a ~ 0.3 as well. (Note that 
in that example as was also implied from the option prices, whereas 
in the current example, as was calculated from historical observations 
for each entity). 

Finally, we show a figure (Figure 1) by way of which we want to 
give a feel for the effect of tails {q) and skew (a) on the credit spread 
values. We used r = 4%, as = 20%, T = 1 year and = 100 
(in arbitrary monetary units). We varied the ratio d = D/Aq from 
0.5 to 1.5 and calculated the spread according to Eq (17) for various 
combinations of q and a, calibrating them by adjusting as so that all 
curves coincide for d = 1. Then, the corresponding volatility which 
would have been needed for the standard Merton model {q = a = I) 
of Eq (9) to reproduce the same values of the spread were backed 
out. These implied volatilities were then plotted out as a function of 
d. Similar in spirit to the implied volatility smiles commonly used to 
depict the deviations that tails and skew have on option prices relative 
to the standard Black-Scholes-Merton model, these curves can give an 
intuition as to how the spread values vary relative to the standard 
log-normal Merton model. 

The standard model corresponds to g = 1 and a = 1. The effect 
of increasing the tails a bit can be seen in the curves corresponding to 
q = 1.2 and a = 1: the standard Merton model undervalues leverage 
ratios greater and less than 1. As q increases (see the curve for q = l.A 
this effect is enhanced. The curve q = 1 and a = depicts the 
effect of skew, and corresponds to a CEV type Merton model. The 
standard Merton model overvalues spreads if there is a debt-to asset 
ratio greater than 1, and undervalues it otherwise. The effect of both 
tails and skew can be seen in the other curves. The observed behaviour 
is consistent with the intuition that fat tails correspond to higher 
volatilities with respect to a Gaussian model, and increased left skew 
would correspond to higher volatilities for d < 1, and relatively lower 
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volatilities for d > 1. 

It is quite easy to visualize where the preliminary empirical results 
of g = 1.3 and a = 0.3 [17] would lie in this curve. However, further 
studies must be done to analyze the CDS values systematically as a 
function of the ratio d. In the initial study, the results reported were 
as average over all the different companies, each one which had its 
own capital structure and leverage ratio. However, encouraged by the 
fact that the parameters which best describe the CDS spreads are in 
the same ballpark as those which well-fit stock options as empirical 
return distributions, it might be worth while to push a little further 
along the path we have proposed, a project we are currently pursuing. 



5 Appendix 

Xf is a function of ilf given by 

and f = (e2(°-i)'-^l)/(2(a - l)r). Pg is given by equation Eq (13) 
evaluated at t = T. The coefficients a, 6, c and d are given by 

a = 5o(g-l) + ^7 (19) 



b = ad — rjgi 
V 

92{q - 1) 



c = (q-l)^d 
V 



d = 



with go = l{f)^^,gi = 7(T)^, 52 = 9^ and ^{f) = ((3 - 

q-l ^ 2 

q){2 — q)cq)^-iT^-'! . The payoff condition At = D yields a quadratic 
equation with the roots 



with 



1 — a 2 
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M = c—-ad (22) 
R - ^ + a— (23) 



These results are as in [8]. 
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